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Abstract 

In this paper we apply the formalism of local composite operators as developed by Ver- 
schelde et al. in combination with a constant chromomagnetic field as considered in the 
seventies by Savvidy and others. We find that a nonzero {A 2 ^) minimizes the vacuum en- 
ergy, as in the case with no chromomagnetic field, and that the chromomagnetic field itself 
is near-to zero. The Nielsen-Olesen instability, caused by the imaginary part in the action, 
also vanishes. We further investigate the effect of an external chromomagnetic field on the 
value of {Ajj), finding that this condensate is destroyed by sufficiently strong fields. The in- 
verse scenario, where (At) is considered as external, results in analogous findings: when this 
condensate is sufficiently large, the induced chromomagnetic field is lowered to a perturbative 
value slightly below the applied {A 2 ^). 

1 Introduction 

In the seventies Savvidy pQ found that the SU(2) Yang-Mills vacuum is unstable against 
formation of a constant chromomagnetic field. This new vacuum, though, is neither gauge 
nor Lorentz invariant. 

Not much later Nielsen and Olesen [2] showed that the action in this new vacuum has an 
imaginary part, meaning that the Savvidy vacuum is unstable as well. Ever since then many 
ways have been explored in order to stabilize this, the most well-known being a dynamical 
Higgs approach [3] E] an d the "spaghetti vacuum" [3] consisting of a superposition of many 
domains with different orientations, forming a kind of liquid crystal. Recently new roads of 
investigation have been explored using the Cho-Faddeev-Niemi decomposition [5j [6] , see for 
example 0HH5]. 

In [10] Gubarev, Stodolsky and Zakharov proposed that the condensate {A 2 ^) might have a 
significance in the Yang-Mills vacuum. Although gauge variant, this quantity has been found 
to be relevant in detecting the condensation of magnetic monopoles in compact QED. It can 
be shown that this quantity is minimized when working in the Landau gauge. Hence, in that 
case, (A^) m i n can be given a gauge invariant interpretation [ID] . 

In [11U12] one of us introduced the formalism of local composite operators (LCOs) so as to 
enable them to calculate this quantity. In [13] their method was applied to SU(N) Yang-Mills 
theory, indeed resulting in a nonzero value for (A^) m i n . A consequence of a nonvanishing value 
for (A^) is the dynamical generation of an effective gluon mass. 

In this paper we will combine the formalism of LCOs with the constant chromomagnetic 
background field of Savvidy. In section [2] we will give a short review of the LCO formalism in 
Yang-Mills theory. Section [3] will be devoted to calculating the effective action, which will be 
discussed in section [4] There we search for minima of the action, and we consider the effect 
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of each field on the induced value of the other one. Finally, in section [S] our conclusions will 
be presented. 



2 LCO formalism 

In this section we will review the LCO formalism as proposed in [13] . 

As a first step the gauge is fixed using the Landau condition, i.e. the linear covariant 
gauge 9 M A M = with £ — * 0. Then, a term 

\z 2 JAl (I) 

is added to the Lagrangian density. Here Z2 is a multiplicative renormalization constant and 
J is the source. As it stands, the theory is not renormalizable. To correct this a new term 

- \Z^J 2 (2) 

has to be added. Here £ is a new coupling constant and Z( is its renormalization factor. This 
Lagrangian is now multiplicatively renormalizable, as shown in [14] using a BRST analysis. 
There are several problems, though. 

As a first problem we have introduced a new parameter, (, creating a problem of unique- 
ness. However, it is possible to choose ( to be a unique meromorphic function of g 2 based on 
the renormalization group equations. In [13] they found using the MS scheme in d = 4 — e 
dimensions (up to one-loop order and with N c the number of colors): 

_ 9 Af 2 - 1 1 N 2 ~ 1 161 
4 ~ 13 N c ff 2 + 16^ 2 52 W 

Secondly the presence of the J 2 term spoils an energy interpretation for the effective 
potential defined via the Legendre transform. In order to solve this, a Hubbard-Stratanovich 
transformation is applied by inserting unity into the path integral: 



1 = Mj [Va] exp-^/g + \Z 2 Al - Z c ( j) 2 *»- 



(4) 



with A/" an irrelevant constant. This eliminates the IZ2JA 2 and Z^J 2 terms from the 
Lagrangian and introduces a new field a. The result is: 

e~ w(J) = /p?yy[2V]exp- / (£ YM [^ M ,c,c] + £ L co[^M,a]--j)d 4 :r. (5) 



y ^£ym[^, c, c] + £lco[^ m , ct] - ^ x 



Herein £ym is the well-known Yang-Mills Lagrangian with Faddeev-Popov ghosts, fixed in 
the Landau gauge, and 

Uoo\A» M = + l^gaAM + \^KK? (6) 

Now J acts as a linear source for the a field, so that we can straightforwardly compute the 
effective action Y(o) using the above expressions. 

If we compare this to the original expression, we find that the expectation value of a 
corresponds to the expectation value of the composite operator 

a = -g(~Z 2 Al-Z < CJ^ (7) 

In the limit J — *■ this operator corresponds (up to a multiplicative factor) to A 2 . We can 
also read off the effective gluon mass in lowest order: 
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3 Effective action 



3.1 Introductory matters 

We now proceed to combine the formalism of a constant chromomagnetic background field 
with the one of massive gluons using the LCO formalism. 

Since we are now working with a background field, it is more appropriate to use the 
Landau background gauge [TS] D^AlA^ = instead of the usual Landau gauge prescription 
d^Afj, = 0. Here is the background field. In order to do so, some alterations are in order. 

A BRST analysis (for BRST in the background gauge, see for example [16]) shows that, 
in order for the LCO formalism to stay renormalizable, the condensate A 2 ^ must be replaced 
by 



At 



An, ~\~ 2 A n A n 



(9) 



with A M the total gauge field and _4 M the quantum fluctuations, A^ — _4 M + A^. 

This replacement will change nothing in the expressions for ( and the renormalization 
constants. In the limit A M = these must reduce to the original expressions, and since these 
constants are dimensionless and A M is the only dimensionful parameter they could otherwise 
depend on, they will not change when switching on a nonzero background field. 

As a result, the action we depart from is given by 



+ 



+ 



, <r/ ^ ^ . 2 , / ^ fC ^(W + 2^) + I^( W 2 + 2^^) 2 (10) 



For simplicity, we will work in SU(2). If we choose the background field S" to be a 
chromomagnetic field in the 2-direction in space and in the 3-direction in isospace, we can 
write 

A; = H Xl 5 a3 g^ (11) 
With this expression, the effective potential at one loop is given by: 



V eS = i^ + ^-logdet^) 

+ ilogdet L^S ab J^ga - g^V 2 ab + (l - ± ) liyP,. i ! - 2< /( -'"' : 7/,s'; 



(12) 



where the limit ^ — + is implied, and 



-1 

1 



(13) 



0/ 



3.2 Spectrum of V 2 

We start by calculating the determinant of the ghost operator. 

If we use the eigenbasis of e abi , we first have the "3" ghosts, for which the covariant 
derivative reduces to an ordinary one, and then we have the "+" and the "— " ghosts with 
eigenvalues ±i. For those last ones the covariant derivative equals D M = <9 M ± igHxig^- 

The "3" ghosts give a trivial contribution of trlog9 2 = 0. 

For the "+" and the "— " ghosts we need the eigenfunctions of T> 2 = d 2 ± 2iHx\gd2 — 
g 2 H 2 x\. This is a harmonic oscillator, and we readily find: 




where ip n (n 6 N) is the nth eigenfunction of the harmonic oscillator and with a = (03,04). 
We get with dimensional (d — 4 — e) and zeta function regularization (}2 n (n + q)~ B — £(s; g), 
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the Hurwitz zeta function): 



Tj +OQ p jd — 2 

trlogD 2 = A^ln^n+D+p 2 ) (15) 



n=0 

9 2 H 



3(4tt) 



0+l-ln^-12C'(-l)-ln(2)) (16) 



where f(s) = X] n n s i s the Riemann zeta function. Here, the factor of two comes from the 
contributions of both "+" and "— " ghosts. 

3.3 Spectrum of the gluons 

The gluons can be split into two classes: the ones obeying the Landau gauge prescription 
£>,^m = 0, giving 

\ tr-p^^o log {s^ ah -^^ 9 a ~ 9iwT%b + 2ge ab3 HSl v ^j (17) 
and the ones not satisfying the prescription, giving 

2 tru M ^ M ^ l°g T>pV v + constant (18) 

where the irrelevant constant part contains lim£_ohi£- 

The spectrum of this last operator can be reduced to the spectrum of D 2 . If i/> M is an 
eigenfunction of T>^T> V with eigenvalue k 7^ 0, we also have that T> 2 D^ip^ = kV^ip^ so that 
©m'/V i s an eigenfunction of T> 2 with eigenvalue k. This means that all eigenvalues of T> ll T> v are 
also eigenvalues of T> 2 . Conversely, if / is an eigenfunction of the operator V 2 with eigenvalue 
p, then 25 M / will be an eigenfunction of T> I1 'D V with the same eigenvalue. Thus we see that 
these two operators have an identical spectrum and we can write 

tr ©^ tI #olog£ , M £ , „ = trlog© 2 (19) 

The expression on the right-hand side has been calculated above. Since the ghosts will come 
with a factor —1 and the gluons with a factor 1/2, exactly minus one half of the result given 
there will remain. 

For the gluons fulfilling the gauge prescription, we start with the "3" gluons. A straight- 
forward calculation yields for the three polarizations in dimensional regularization: 

e ,__ ( Z2 2\ 3Z 2 a 2 (2 5 , Z20 



2 ■ w n« F "t"wFt + r h w l (20) 



Secondly there are the "+" and the "— " gluons: 
1 , ( Z 2 



| tr-D t ^ ti=0 log I gM^ p^ g ff - g^v 2 ± iigHSl u J (21) 

with D M = <9 M ± igHxig^- We start with the "+" gluons. We now pass to the polarization 
basis wherein St u is diagonal. We get: 



o3 _ 



—i 



V 0/ 



£V = {VgHa iyfg~Ha< d 3 9 4 ) (22) 



with a and the ladder operators of the harmonic oscillator from equations (|14|) . This reduces 
the problem to four one- dimensional harmonic oscillators with the same eigenfunctions as in 
([Till. Tne eigenvalues are gH(2n + 1 + 2s) + p§ + pi + Z 2 cr/{gZ c C;) with s = -1, 1, 0, the 
spin eigenvalue of the state. 
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Now we have to restrict the spectrum according to the Landau background gauge. For 
this purpose we construct the following vector functions from the scalar eigenfunctions 



fn = e 




gHxi 

3** + vm) J 



_P2_ | 

sfgH) 



n = -1,0, 1,2.. 



(23) 



where tp n with n negative is defined to be zero. The vector e M is a polarization vector. These 
functions have eigenvalues Z 2 a j {gZ^C,) +gH(2n+ 1) +p§ + p\. To see whether they obey the 
gauge condition, we calculate 



n = -\ 

(ei V / p77 + ie 3 p 3 + ie 4 p 4 ) e xS ' p tp (. . .) n = 

_ (e 1 ^/g~H^n + 1 + ie 2 *fgl[*~fn + ie 3 p z + ie 4 p 4 ) e la: ' p i/>„(. . .) n > 



(24) 



We conclude that, for n = — 1, there is but one polarization with a contribution of 1/2 tr \og(Z 2 (J / {gZ^ ()- 
gH + p 3 +P4). For n = 0, of the three, one is eliminated by the gauge prescription, leaving us 
with 2 polarizations (2— e in dimensional regularization) each contributing 1/2 tr log(Z 2 a / (gZ^Q + 
gH + P3+P4). For n > we have the usual 3 (3 — e) polarizations with the usual contribution. 
For ease of calculation, we calculate the second and third groups together with 3 polarizations, 
so that we have to subtract the contribution of n = exactly once. 
The gluons with n = — 1 give 



gH f d 2 - e p 
4^ J (2tt) 2 - £ 



log 



\gz<<; 

For the gluons with n = there remains: 



gH 

47T 



! — -^ + pJ= ^ | v 7 + 1 - ln /i 2 j (25) 



(4tt) 2 



+ 1 - In 



+ gH 



M 2 



And finally all the other states contribute 



(26) 



9 2 ZfC 2 



g 2 H 2 



4(4tt) 2 



2 ffj fA 6(gH) 2 d( 



3 ln / + (4tt) 2 



9s 



1 

' 2 2g' 2 HZ c C 



(27) 



Here, </) denotes the analytic continuation of the Hurwitz zeta function, which for first 
argument greater than one is defined as 



fc=i 



or by its integral representation 



1 f°° t"-\ 



-dt 



(28) 



(29) 



The derivative dt^/ds stands for the derivative with respect to the first argument. In this 
last calculation we have made use of the relation between the Hurwitz zeta function and the 
Bernoulli polynomials, in our case: 

C(-l,*) = -^ = -£ + £-i (30) 



The "— " gluons give exactly the same contribution, so that the above expressions must 
be multiplied by a factor of two. 
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Figure 1: The effective action (|31[1 in function of H and a' . The black line denotes a' = H, where 
the imaginary part of the action vanishes. 



3.4 Total 

If we sum all the terms we have calculated, and we substitute the values for the renormalization 
constants, we get: 

v - V + 27a ' 2 

9 g V 2 / 1 161 1 9^__2 2gH \ 
~4(4tt) 2 ^2 + 78 3 % 2 3 n p? ) 
2g 2 Ha' a' - H 
(4^p~ n a' + H 

-^(4-21n^-21n^ + Iln^-2CVl)-iln(2) 

+ (4tt)2 9s ^ i; 2 + 2£fy! 
where we have set 

9(7 = 27 5 °" ^ ' 

so that the effective gluon mass squared is m 2 ff = ga' . The real part of (|31[1 is plotted in 
figure [1] 

In the limit — *• this expression reduces to the one obtained in [13], and when taking 
o — > we get the result of Nielsen and Olesen [2] modulo some differences due to the use of 
another subtraction scheme and gauge. 

When H > a' our potential (|31l) has an imaginary part 

(33) 

which reduces to the Nielsen and Olesen result for a' = 0. It turns out to vanish for H < a' , so 
that the Nielsen-Olesen problem of the Savvidy vacuum is then resolved, as predicted would 
happen by Nielsen and Olesen [2] themselves. 

4 Discussion 

In the next two subsections, we will find that the minimum of the effective potential is for 
H — (or virtually zero) and a' the value calculated in [13] . 
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Figure 2: The difference between E va , c (H) and E V!!IC (H = 0) for very small values of H in the non- 
perturbative minimum for a' . A shallow minimum (order 1(T 27 Afjg-) is seen for H » 7 x KT 13 AHjg. 



In order to do so, we will consider two cases: first H will be considered as an external field 
and a 1 as an effective gluon mass induced by quantum effects, and next we will investigate the 
influence of a nonzero a' on the value of the Savvidy field. When looking at small values of the 
external fields, the analyses can be done analytically by expanding the potential in this small 
parameter. The scale can then be chosen according to renormalization group considerations. 
For bigger values of the external fields, however, we proceed numerically. In this last case 
the scale p? is, for the ease of calculation, fixed equal to p? — 4.12A^g-, the value of ga' in 
the global minimum of the effective action. In that point the coupling constant is reasonably 
small: 

g 2 36 



8.3 " 187 -°- 19 (34) 



4.1 Effect of H on a 1 



If H is set to zero, the effective potential has a perturbative extremum (a maximum) in a' p — 
and a non-perturbative minimum at 

S <, = Afjge^ =4.12A| i g (35) 

where the scale was chosen equal to ga' np . 

For small H the equations can be expanded in a series in H: 

v m '\ 27 a ' 2 9 ff V2 f 5 i 161 , 9°' 
V e a{H,a ) = jg-g 4(4^)2 + ~~yg ~p~? 

+ 0{H 3 \tiH) (36) 

To obtain this, we used the expansion of the Hurwitz zeta function for large arguments, which 
can be straightforwardly calculated from the integral representation (|29|) . From this can easily 
be obtained that, up to this order, 

ga' np = gcr' H=0 - — (37) 

ycr H=0 

so that a' np decreases with a raising of H . The vacuum energy changes like 

*- = E »- - h" 2 + w? ( w + \ + 2 ^<- 1) ) + (38) 

We see that for very small H the term of order H 2 In H will dominate, lowering the vacuum 
energy. Very fast, though, this term will be supplanted by the terms of order H 2 , and the 
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Figure 3: Real (full line) and imaginary part (dashed line) of the potential for H = 0.4 A^-. 



energy will start increasing again. The effective potential in this regime is depicted in Figure 
[2] The lowest value is reached when 

H = a'exp(^^- -lgZ_21n2-12C'(-l)') = 3.92 x 1CT 13 a' (39) 
V 13<7 2o / 

Since this result is astronomically small, there is no reason why it wouldn't disappear when 
higher-order corrections or any other effects are taken into account. For all practical purposes 
one can say that the vacuum energy is lowest when H = and a' has the value given in (I35|) . 
One would expect terms containing In// to be replaced with ln(// + a) when switching on a 
mass, causing this residual chromomagnetic field to vanish. This does not happen, though, 
because the ghosts and the unphysical gluon do not cancel. This is related to the unitarity 
problem of the model, which could be solved non-perturbatively in the zero color sector when 
incorporating confinement. [13] 

When H is increased, analytic methods have to be abandoned and we solve the equations 
numerically instead. A qualitative sketch of the effective potential in this regime is depicted 
in Figure [3] A nonzero imaginary part exists ever when a' < H, as mentioned above. In 
the real part of the action, the point with a' = is no longer an extremum, but a new 
perturbative minimum forms for a' between zero and H. This is separated from the original 
non-perturbative minimum by a little hill with a top at a' slightly above H . The value of 
a' in the non-perturbative minimum decreases with increasing H. For higher H a point is 
reached where the minimum with smaller a' has a lower energy than the one with greater 
a'. We thus find a first-order phase transition around H = 0.40 A^g-. For H yet higher, 
the non-perturbative minimum disappears altogether and only the perturbative one remains. 
These evolutions can be seen in Figure [4] 

The conclusion is that a nonzero chromomagnetic field decreases the effective gluon mass, 
and when the field is sufficiently high a phase transition occurs, lowering the mass to a value 
slightly lower than gH. 

4.2 Effect of a' on H 

We can take the limit a' = 0, giving 

Here we used that £(s, 1/2) = (2 s — l)f (s). Ignoring the imaginary part, and putting p 2 equal 
to the value of gH in the global minimum, we obtain a perturbative extremum in H = and 
a non-perturbative one in 

9H = p 2 exp (-gj + g + |f C'(-l) + ln(2)) - 1.71 A^ (41) 
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Figure 4: Left: The various values of a' as functions of H . The full line is the non-perturbative 
value of cr', the dashed line is the value of cr' in the lower minimum, and the dotted line is a' — H , 
drawn for reference. Right: The vacuum energy in the minima as a function of H. In both plots 
the thick dot indicates where the higher minimum in the effective potential disappears. 




Figure 5: Real (full line) and imaginary part (dashed line) of the potential for a' — 0.5 A| 



When expanding in a series in cr', the next term in the real part is 

3g 2 Ha' 



(4tt) 2 



•In 2 



(42) 



meaning the non-perturbative minimum will be lowered. This also increases the induced value 
of H by an amount of 9a'ln2/22 « 0.28 cr'. 

When going to higher values of cr', we find that H — (or near-to zero) turns into a local 
minimum of the potential. For H slightly below a' there is a maximum and for H higher than 
cr' there is a non-perturbative minimum (see Figure [5}. When increasing cr', this last one first 
deepens out, reaching a lowest value for a 1 = 0.40 Afjg-, and it then goes up again. The value 
of H in this point grows with increasing cr'. For cr' big enough this H asymptotically goes to 
a'. 

The value of the effective action in H = decreases for rising a', so that around a' — 
0.48 A|jg- it dives lower than the energy in the non-perturbative minimum. This means that, 
at this point, there is a first-order phase transition from the state with H > a' to the one 
with H « 0, causing the imaginary part in the action to vanish. This is depicted in Figure [6] 

We conclude that switching on a nonzero gluon mass first makes H increase, and then 
destroys it completely. When the gluon mass is sufficiently large, the vacuum is no longer 
unstable against the formation of a constant chromomagnetic field, and the Nielsen-Olesen 
instability, caused by the imaginary part, also is resolved. 
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Figure 6: Left: The induced value of H as a function of a' (full line). For higher values of a' this 
nears the asymptotic H — a' (dashed line). Right: The vacuum energies in the non-perturbative 
minimum (full line) and in H — (dashed line). The branch for H = is the same as the potential 
calculated in [13], reaching its lowest value for a' np = 1.06 A™-. 

5 Conclusions 

We found that, when considering both a constant chromomagnetic field and an (A 2 ^) conden- 
sate, the effective action was minimized for zero (or near-to zero) chromomagnetic field with 
a non-perturbative value for (j4?) as found by Verschelde et al. [13] There are no unstable 
modes any longer, and the imaginary part in the action is zero in this minimum. 

When considering the situation in which H is an external field, we found that applying 
such a field first lowers the value of the induced mass, and for H around 0.40 Ajt^ the non- 
perturbative mass is destroyed, leaving only a perturbative value slightly smaller than H. 
The action then has a (small) imaginary part as in the Savvidy case. 

When, on the other hand, considering the effect of the mass on the Savvidy field, we found 
that a sufficiently high gluon mass (a > 0.48 A|jg) destroys the induced H field, at the same 
time causing the Nielsen-Olesen instability (the imaginary part in the action) to vanish. 
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